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ABSTRACT 

In the past decade quantum algorithms have been found which outperform the best classical 
solutions known for certain classical problems as well as the best classical methods known 
for simulation of certain quantum systems. This suggests that they may also speed up the 
simulation of some classical systems. I describe one class of discrete quantum algorithms 
which do so — quantum lattice gas automata — and show how to implement them efficiently 
on standard quantum computers. 
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1. Introduction 

Quantum computing originated with Feynman's observation that quantum systems are 
hard to simulate on classical computers, but that they might be easier to simulate if one 
had a computer which could be operated quantum mechanically [1]. Developments dur- 
ing the subsequent two decades have not only supported this observation [2-8], but have 
also demonstrated that quantum computers would — if they existed — solve certain classical 
problems like factoring more efficiently than is possible on classical computers running the 
best classical algorithms known [9]. This raises a natural question [10-12]: Might quan- 
tum computers help with the simulation of classical systems? Or more specifically, given 
the focus of this workshop, could quantum computers efficiently simulate fluid dynamics? 
Without answering the specific question, in this paper I try to explain why the answer 
to the general question may be 'y^s', and along the way explain some of the quantum 
algorithmic tricks which seem likely to be useful in future investigations of these questions. 

I begin by (very) rapidly introducing the ideas of quantum computation in §2. For 
a complete presentation, see [13]. In §3 I describe one approach to simulating quantum 
systems — quantum lattice gas automata (QLGA) [14] — and then in §4 explain a connection 
with simulation of classical systems. The crucial issue is the relative complexity of quan- 
tum and classical algorithms; so §5 contains a detailed analysis for one specific problem, 
including some new results on implementing QLGA algorithms on 'standard' quantum 
computers. I conclude with a brief discussion in §6. 

2. Quantum computers 

The possible states of a classical computer are (very long) bit strings h\ . . .hn G {0, 1}"^. A 
particular computation proceeds via a sequence of maps to new bit strings: h'^. . . b'^, .... A 
fundamental result, which contributed directly to the conceptual development of quantum 
computation, is that any classical computation can be made to be reversible [15], i.e., these 
maps can be chosen to be permutations on the space of states. Quantum computation can 
then be understood as a generalization of classical computation: The possible states of a 
quantum computer are superpositions of bit strings ^ a6i...6„ \bi . . . hn) G (C^)®" (each 
tensor factor is called a quantum bit or qubit), where ^ |o6i...6„|^ = 1 so that the norm- 
squared of each amplitude ai,^...i,^ e C is the probability that the state |6i . . . 6„) is observed 
if the quantum system is measured in this 'computational basis'. A particular quantum 
computation proceeds via a sequence of unitary maps to new states X^'^fei 6„l^i • • -^n), 
. . .. This much is a generalization of classical reversible computation since permutations 
are unitary maps, and each classical state is an allowed quantum state. The difference 
is that the final state is not directly available; it can only be sampled according to the 
probabilities given by the norm-squared of the amplitudes. 

To evaluate the computational complexity of an algorithm, either classical or quan- 
tum, we must specify a set of elementary operations, the number of which used during 
the computation quantifies the complexity. If we allow arbitrary permutations classically, 
or arbitrary unitary transformations quantum mechanically, any state can be reached from 
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any other state in a single step — these are 
clearly not reasonable models of compu- 
tation. Instead, consider the 'gate opera- 
tions' shown in Fig. 1: NOT, C-NOT and C- 
C-NOT ('C' abbreviates 'CONTROLLED'). 
Each of these is a permutation on the set 
of bit strings, and Toffoli [16] has shown 
that this set of three gates is universal for 
classical computation, in the sense that 
any (reversible) boolean operation can be 
decomposed as a sequence of these gates — 
which by Bennett's result [15] suffices for 
universality. Each of these gates can be 
extended by linearity to a unitary map on 
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acting non-trivially only on a sub- 



Figure 1. A universal set of classical reversible 
gates. Inputs are on the left; outputs are on the 
right. + denotes addition mod 2. 



set of 1, 2, or 3 qubits, and thus can also 
be used as a quantum gate. (Tradition- 
ally the quantum NOT gate is denoted by 

ax or X.) Two other quantum gates which are particularly useful act on single qubits: 



"^7=2il -1 



and Ru 



1 
CO 



the 'Hadamard transform' and phase u rotation, respectively. These matrices have been 
expressed in the computational basis; thus 

iy|o)^-L(|o) + |i)) 
^|i)-^(|o)-|i)). 

C-NOT and Rcj for uj = e*^ with 9 = cos"-*^ | form a universal set of gates for quantum 
computation [17]. 

As we noted earlier, any array of (reversible) classical gates can be simulated by some 
array of quantum gates. The remarkable fact is that in some cases fewer gates are required 
quantum mechanically. The following example is the smallest version of the Deutsch-Jozsa 
[18] and Simon [19] problems: 

Example. Given a function / : {0, 1} {0, 1}, we would hke to evaluate /(O) XOR 
/(I). The function is accessed by calls which have the effect of taking classical states (x, b) 
to (x, b © f{x)) where © denotes addition mod 2. This is a reversible operation and thus 
also defines a unitary transformation on a pair of qubits: /-C-NOT. Classically, this gate 
must be applied at least twice (once with a; = and once with x = 1) in any algorithm 
which outputs /(O) XOR / (I) correctly with probability greater than ^ (assuming a uniform 
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distribution on the possible functions). Quantum mechanically we can exploit interference 
to do better, applying the operation only once. Suppose the system is initialized in the 
state |0) |0). Then apply the following sequence of unitary operations: 

|O)0|O)^^ ^i|a;)®(|0)-|l)) 

x=0 
a;=0 

E ^(-i)^'^^n-i)^iy)®(io)-ii)) 

x,y=0 ^ 

= i=|/(0) XOR/(1))®(|0)-|1)). 

The first step Fourier transforms the 'query' qubit into an equal superposition of |0) and 
|1), and initializes the 'response' qubit into a state which will create the phase (— 1)-^'^'^) 
in the second step. The third step Fourier transforms the query qubit again, creating the 
interference which ensures that subsequently measuring the query qubit outputs /(O) XOR 
/(I) correctly with probability 1. The gate array implementing this quantum algorithm is 
shown in Fig. 2; it includes only one /-C-NOT gate. 

|0> [e] [T\ [h} If (0) XOR f (i)> 

|0> © I0> - |i> 



Figure 2. The gate array implementing the quantum algorithm which solves the XOR problem 
with a single call to the /-c-not gate. The states are not normalized. 



This example displays perhaps the simplest quantum improvement over classical com- 
putation. Other superior quantum solutions to classical problems include Deutsch and 
Jozsa's balanced function algorithm [18], Bernstein and Vazirani's parity search algorithm 
[20], Simon's period finding algorithm [19], Shor's factoring algorithm [9], Grover's unstruc- 
tured search algorithm [21], van Dam's algorithms for weighing matrices and quadratic 
residues [22] , and Hunziker and Meyer's highly structured quantum search algorithms [23] . 

3. Quantum lattice gas automata 

Although the example in §2 demonstrates the superior computational power of quantum 
systems for certain problems, it seems to have little to do with confirming Feynman's 
original idea that quantum systems could efficiently simulate other quantum systems [1]. 
For the community attending this workshop, a natural place to look for such confirmation 
is lattice models: The possible configurations for each particle on a one dimensional lattice 
L are labelled by pairs {x, a) E L x {±1}, where x is the position and a the velocity. 
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A classical lattice gas evolution rule consists of an advection stage (x, a) t— > (x + a, a), 
followed by a scattering stage. Each particle in a quantum lattice gas automaton (QLGA) 
[14] exists in states which are superpositions of the classical states: = ^ 'i/'x,a|2;, «), 
where 1 = (V'lV') = ^ V'x.aV'x.a- The evolution rule must be unitary; the most general 
with the same form as the classical rule is: 

^ Vx,a|a:, a) h^:^^^ + a, a) 

scatter \ i ri \ , /\ 

where the scattering matrix is 

g _ f cos s i sin s \ 

y i sin s cos s J 



Let U denote the complete single timestep 
evolution operator for a single particle, 
the composition of advection and scatter- 
ing. Fig. 3 illustrates this quantum evo- 
lution: at s = it specializes to the clas- 
sical deterministic lattice gas rule. The 
Ax = At limit of this discrete time 
evolution is the Dirac equation [14]; the Ax^ 
[24]. 




Figure 3. The general evolution rule for a single 
particle in the one dimensional QLGA. 

= At — s> limit is the Schrodinger equation 



QLGA models can be generalized to higher dimensions [24] , and to include more parti- 
cles [14,25,24], potentials [26,27] and various boundary conditions [28]. These are quantum 
models which we might try to simulate classically or quantum mechanically. Figures 4-6 
show the results of classical simulations of plane waves, wave packets, and scattering off po- 
tential steps, respectively. These support the claim in the previous paragraph that QLGA 
are discrete models for quantum particles. In the next section we will see how they are also 
relevant to the question of simulating classical physical systems quantum mechanically. 



Re 




Re 



Figure 4. Plane waves in the general one dimensional QLGA [26]. 



5 



Quantum computing classical physics 



David A. Meyer 




Figure 6. Scattering of a wave packet from a potential step in the general one dimensional 
QLGA [26]. 
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4. Diffusion 

The evolution rule for a single particle QLGA bears some resemblence to a random walk. 
More precisely, it is the unitary version of a correlated random walk [29,30] — with = tt/A 
the analogue of uncorrelated. In Fig. 5, for example, we can see evolution like that of a 
biased random walk, with a spreading Gaussian distribution. And, in fact, this is as close 
as possible — there is no quantum random walk in the sense of local, solely x dependent 
unitary evolution; the a dependence must be included [31]. (This is not true for quantum 
processes continuous in time [32] , which can also be used for computational purposes like 
the ones we are considering here [33,34].) 

Consequently, there are differences. Diffusion approaches an equilibrium state, in- 
dependently of the initial condition (on connected spaces). Unitary evolution of a single 
particle QLGA cannot: the distance 5 between successive states, defined by cos 5 = {Uiplip) 
is constant, so the evolution of any state which is not an eigenstate of U with eigenvalue 
1 does not converge. Each state implies a probability distribution on the lattice, given by 

Pt{x)=proh{x;t) = |^^,_i(t)|2 + |V^,,+i(t)|2, 

so we can also ask if this converges. In fact, this probability distribution is constant for 
each eigenstate of U [26]. But since there exists T G Z>o such that Af ~ 1 for all 
eigenvalues Aj of U, and hence U'^ip ~ ip, for any initial state such that Pi ^ Pq, the 
probability distribution cannot converge either [35] . 

Aharonov, Ambainis, Kempe and Vazirani have shown, however, that the time average 
of the probability distribution does converge [35]: Expand the initial state = 
in terms of the eigenvectors of U . Then U^\i\)) = ajA*|0i), so 

I V 2 

pmh{x,a;t) = |^ ajA-(a;, Q;|0i) 

i 

= ^ aiaj{XiXjY{x, a\(l)i){(l)j\x, a). 

Then the time average of the probability is 

^ T-l ^ T-l 

-^pToh{x,a;t) ^^aiaj{x,a\(l)i){(j)j\x,a)- ^{XiXjY. 

t=0 i,j t=0 

For Xi 7^ Aj, the interior sum goes to as T ^ oo. This leaves only the terms in the 
sum for which Xi = Aj, which are independent of T. Thus the time average converges. 
In particular, for the one dimensional single particle QLGA, it converges to the uniform 
distribution which is the equilibrium distribution for diffusion on one dimensional lattices. 
That is, by measuring the position at random times we can simulate sampling from the 
equilibrium distribution of classical diffusion. Although this is not true for all graphs. 
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e.g., the Cay ley graph of the nonabelian group S3 [35], we have analyzed one example of 
discrete quantum simulation of a classical physical process. 

5. Computational complexity 

To be truely useful, the quantum computation should be more efficient than the corre- 
sponding classical computation. Classically, 0{N'^) steps of a random walk are required 
to approximate the equilibrium distribution on a lattice of size A^. Aharonov et al. have 
shown that only O(A^logA^) steps of the single particle QLGA are required for equally 
accurate sampling [35], and more detailed calculations by Ambainis, Bach, Nayak, Vish- 
wanath and Watrous [36] show that 0{N) steps suffice. The proofs of these results de- 
pend on careful estimates about, for example, the distribution of eigenvalues of U and 
are somewhat involved. Simple simulations, however, provide a heuristic explanation for 
this quantum improvement. Fig. 7 shows the evolution of a QLGA particle initialized 
at — (|0, — 1) + |0, +l))/\/2. Notice that the peaks of the probability distribution — 
indicated by the darkest squares in the plot — spread approximately linearly in time. This 
is the origin of 0{N) number of steps required for the QLGA to sample the equilibrium 
distribution. 

Thus the QLGA provides a quadratic improvement in the number of steps required 
compared to a classical random walk. To verify a computational improvement, however, 
the QLGA must not require much more computation per step. Each step of the random 
walk requires a coin-flip, i.e., generation of a random number. Inasmuch as this is possible 
classically, it requires constant time, independent of N. Each step also requires the addition 
of the ±1 result of the coin flip to the current position. Since the latter is a log N bit integer, 
this requires 0(logA'") elementary operations. Thus the total number of computational 
steps to simulate one run of the random walk is O(A^^logA^). We can compare this with 
the computation required to compute the whole probability distribution by evaluating the 
Markov process. Although this could be computed by matrix multiplication at each step, 
the locality of the process means that we need only compute 

prob(a;; t+ 1) = -prob(a; — l;t) + -prob(a; + 1; t) 

for each lattice point. Since this has a constant computational cost per lattice point, 
evaluating the whole probability distribution for 0{N^) steps takes 0{N^) elementary 
operations. 

For the QLGA we cannot run single trajectories since that would miss any interference 
between trajectories. Thus classical simulation of the QLGA must be like the Markov 
process calculation of the whole probability distribution of the random walk. Again, by 
locality, each step requires constant computation per lattice point. Thus evolution of the 
whole state for 0{N) steps takes 0{N'^) elementary operations. Taking the time average to 
approximate the diffusive equilibrium distribution requires another 0{N) factor, hence the 
same 0{N^) elementary operations we found in the previous paragraph. Unsurprisingly, 
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Figure 7. Evolution of a single particle QLGA initialized at a; = in an equal superposition 

of left and right moving states [14]. Each square corresponds to one of the classical states; thus 
each lattice site is represented by two squares, one for each allowed velocity. The boundary 
conditions are periodic and time runs upward. 
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therefore, we have not discovered a faster classical algorithm for simulating the equilibrium 
distribution of diffusion. 



Thus, to realize a quantum mechani- 
cal improvement we must be able to com- 
pute a single step U of the QLGA evo- 
lution with fewer than 0{N) elementary 
operations on a quantum computer. 
Schematically, U is implemented as shown 
in Fig. 8, acting on log A?" qubits which 
encode the position on the lattice and a 
single qubit which encodes the velocity. 
The scattering operation S acts on the 
last qubit and thus has constant cost per step, 
tiplication by a matrix which is diagonalized b 
left shift \x) I— > |a; — 1) is 





ADVECT 






















Figure 8. The high level circuit for U. There 
are log qubits labelling the position, and a last 
qubit encoding the velocity. 

Advection is a shift operation, i.e., mul- 
'• the discrete Fourier transform F^- The 
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Here oo = e^'^'^l^ and the right shift diagonalizes to the same matrix using rather than 
-Fjv- Assuming A = 2"^, the diagonal matrix can be implemented by n = log A single qubit 
phase rotations: R^^ ® R^^i ® R^^a. ® ■ ■ ■ ® R^^n/2. 

Naively, the Fourier transform Fn re- 
quires 0{N'^) operations. The classical 
FFT improves this to O(AlogA). And 
perhaps the most fundamental result in 
quantum computation is that Fn can be 
implemented with O(log^A) elementary 
quantum operations [9,37]. Thus we can 
implement the advection step with an a- 
controUed-shift operation \x,a) h- > \x + 
a, a), where the shift (in each direction) 
is decomposed as in (1). Slightly less effi- 
ciently, but with the same 0(log^ A) com- 
plexity, we can implement it by acting by a left shift, and then acting by a two step right 
shift, conditioned on ct = 1; this circuit is shown in Fig. 9. Thus on a quantum computer 
we can sample the diffusion equilibrium with 0(Alog^ A) elementary operations. This is 
an improvement by O (A/ log A) over the classical random walk algorithm. 















-F+ 










F - 













Figure 9. A low level implementation of U. The 
phase rotations are as in the text, namely R^^ (8) 
R^2 (8)-Ru;4 each of the three timesteps at which 
they are applied. 



10 



Quantum computing classical physics 



David A. Meyer 



6. Discussion 



We have seen that a classical physics problem — sampling from the equilibrium distribution 
of a diffusion process — is solvable more efficiently with a quantum computer than with 
a classical one, in the sense that the QLGA algorithm outperforms the random walk 
algorithm. The solution utilizes two of the fundamental quantum speedups: First, the 
quadratic improvement in the number of steps necessary for a QLGA simulation compared 
to a random walk simulation is reminiscent of the quadratic improvement of Grover's 
quantum search algorithm — which he, in fact, describes as diffusion [21]. Second, the 
exponential improvement of the quantum Fourier transform over the FFT [9,37] provides 
a speedup in the advection step. Perhaps the simple problem considered here will inspire 
application of these techniques — or new ones — to speed up computation of other classical 
systems with quantum algorithms. For example, a natural generalization would be to 
analyze diffusion in a potential, for which a QLGA simulation has already been shown 
to evolve the mean of the distribution quadratically faster than does a classical (biased) 
random walk, in certain cases [30]. Successful development of such quantum algorithms 
should provide additional impetus to efforts towards building a quantum computer [38]. 
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